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Abstract: Groundwater is a vital water resource in arid and semi-arid areas. Diurnal groundwater table 
fluctuations are widely used to quantify rainfall recharge and groundwater evapotranspiration (ET,). To 
assess groundwater resources for sustainable use, we estimated groundwater recharge and ET, using the 
diurnal water table fluctuations at three sites along a section with different depths to water table (DWT) 
within a wetland of the Mukai Lake in the Ordos Plateau, Northwest China. The water table level was 
monitored at an hourly resolution using a Keller DCX-22A data logger that measured both the total 
pressure and barometric pressure, so that the effect of barometric pressure could be removed. At this study 
site, a rapid water table response to rainfall was observed in two shallow wells (i.e. Obs1 and Obs2), at 
which diurnal water table fluctuations were also observed over the study period during rainless days, 
indicating that the main factors influencing water table variation are rainfall and ET. However, at the 
deep-water table site (Obs3), the groundwater level only reacted to the heaviest rainfalls and showed no 
diurnal variations. Groundwater recharge and ET, were quantified for the entire hydrological year (June 
2017—June 2018) using the water table fluctuation method and the Loheide method, respectively, with 
depth-dependent specific yields. The results show that the total annual groundwater recharge was 
approximately 207 mm, accounting for 52% of rainfall at Obs1, while groundwater recharge was 
approximately 250 and 21 mm at Obs2 and Obs3, accounting for 63% and 5% of rainfall, respectively. In 
addition, the rates of groundwater recharge were mainly determined by rainfall intensity and DWT. The 
daily mean ET, at Obs1 and Obs2 over the study period was 4.3 and 2.5 mm, respectively, and the main 
determining factors were DWT and net radiation. 
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1 Introduction 


In (semi-) arid regions, groundwater plays a significant role in meeting urban, industrial, and 
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agricultural water requirements because of the unreliability or under-allocation of surface water. 
Globally, more than two billion people depend on groundwater, and approximately 50% of 
irrigation water is obtained from groundwater (Alley et al., 2002; Siebert et al., 2010). Due to 
unsustainable development practices, groundwater depletion has taken place in many regional 
aquifers worldwide, such as the North China Plain (Feng et al., 2013), the US High Plains 
(Scanlon et al., 2012), and the Nubian sandstone aquifer in North Africa (Doll et al., 2014). 

Groundwater budget analysis is a critical component for the sustainable use of groundwater 
resources (Bredehoeft, 1997; Zhou, 2009) and the main components of groundwater budgets are 
rainfall recharge and groundwater evapotranspiration (ET,) (Crosbie et al., 2005; Zhang et al., 
2016). Many methods are available for estimating the infiltration recharge from rainfall and ETg, 
including soil water balance (Allison et al., 1994; Gribovszki, 2018), water table fluctuations 
(WTFs) (Cheng et al., 2017; Delottier et al., 2018; Wang et al., 2019a) and numerical simulations 
(Loheide et al., 2005; Liu et al., 2014). Among these methods, the WTF (water table fluctuation) 
method is considered to be effective for quantifying rainfall recharge and ET, because of its 
simplicity and ease of use (Healy and Cook, 2002). The application of WTF to estimate 
groundwater recharge can be traced back to the 1920s (Meinzer, 1923). The procedure used to 
estimate daily ET, using diurnal WTFs was first developed by White (1932), but this method has 
only been widely used in recent decades because of the availability of automatic high-resolution 
groundwater data loggers that can detect millimeter-scale variations (Wang et al., 2019b; 
Gribovszki et al., 2010; Diouf et al., 2020). With recently applied modifications, the sub-daily 
ET, can be estimated from diurnal WTFs (Gribovszki et al., 2008; Loheide, 2008; Yin et al., 
2013). 

Water tables in wetlands are shallow because of the convergence of surface runoff and lateral 
groundwater flows (Toth, 1962; Fan et al., 2013). Consequently, plant species are abundant 
around wetlands because of the high availability of water for their use (Busch et al., 1992; 
Sargeant and Singer, 2016). However, unsustainable groundwater use may degrade wetland 
ecosystems. For example, groundwater-dependent willows died after two years of groundwater 
pumping in a wetland in the USA (Scott et al., 2000). 

The Ordos Plateau is located in Northwest China, where approximately 400 lake wetlands are 
present (Yin et al., 2011a). Natural resources are extremely abundant in the Ordos Plateau, 
including petroleum, natural gas, coal, and sandstone type uranium deposits (Chang et al., 2006). 
It contains the greatest natural gas reserves and the second greatest coal reserves in China (Dai et 
al., 2006). In recent years, groundwater from the area has been pumped out in increasing volumes 
to sustain the rapid development of the Ordos National Energy Base (Hou et al., 2008). The 
effects of overexploitation have already been observed in some areas, resulting in the drying up of 
artesian wells, degradation of groundwater quality, and dieback of groundwater-dependent 
ecosystems (Liu et al., 2015). This is due to the improper quantification of water budget 
components that adversely affect the sustainable development of groundwater resources and 
groundwater-dependent plants in lake wetlands (Liu et al., 2017). In previous studies on the 
groundwater budget in the Ordos Plateau (Wang et al., 2010; Liu et al., 2018), groundwater 
recharge was determined to be the product of annual rainfall and the coefficient of rainfall, 
without considering the impact of rainfall intensity, antecedent rainfall, and the depth to water 
table (DWT) on groundwater recharge. ET, was calculated based on the linear relationship 
between the evapotranspiration (ET) flux at the surface and DWT. However, numerous studies 
indicate that this relationship is more complex than what can be effectively represented by a 
linear model (Baird and Maddock, 2005; Doble and Crosbie, 2017). 

To improve water budget estimates, it is vital to quantify rainfall infiltration and groundwater 
evapotranspiration, as well as their main controlling factors. Therefore, hourly water table 
variations were monitored at three sites over a hydrological year in the semi-arid wetland around 
a lake in the Ordos Plateau. We estimated groundwater recharge and ET, based on the monitoring 
data, and analyzed their controlling factors, such as DWT, rainfall intensity, and antecedent 
rainfall. The objectives of this study were to: (1) explore the response of the water table to rainfall 
and ET, at different depths along a section; (2) estimate groundwater recharge and ET, using a 
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depth-dependent specific yield; and (3) determine the hourly, daily, and monthly variations of 
rainfall infiltration and ET ,, as well as to quantify their controlling factors to explore the 
mechanisms affecting rainfall infiltration and ET, at different time scales. 


2 Materials and methods 


2.1 Study site and experimental setup 


The study site (39°17'26"N, 108°49'17”E) was located within the wetland of Mukai Lake in the 
Ordos Plateau, Northwest China (Fig. 1a). Geographically, it is located in the Middle East Ordos 
Plateau, where the climate is semi-arid and characterized by low rainfall and high ET,. The study 
site was located within an area where groundwater recharges lakes (Yin et al., 2011a; Liu et al., 
2015). The phreatophyte near the lake was native grasses, mainly Carex stenophylla and 
Puccinellia tenuifolia, with a root depth of approximately 80 cm. Although groundwater had not 
been extracted from this lake wetland, an adjacent lake wetland area was heavily pumped, and 
corresponding adverse environmental effects have been observed there, such as wetland shrinkage 
and ecosystem degradation (Liu et al., 2015, 2018). The local government is planning to reassess 
safe pumping rates and requires a detailed water budget analysis. Because of the similar 
hydrogeological and meteorological conditions, the results of this study can provide a reference 
for other lake wetlands in the Ordos Plateau. 

Piezometers were placed in holes drilled with a hand auger at three observation sites, Obs1, 
Obs2, and Obs3, to form a transect perpendicular to the lake (Fig. 1b). The land surface slope and 
hydraulic gradient along the section were 1.7% and 1.0%, respectively. The depth of the 
piezometers varied between 2 and 4 m, depending on the DWT of the site. Each was constructed 
using a PVC pipe with a diameter of 50 mm, screened over the entire length of the subsurface. 
Augured sand was backfilled around the piezometers, and granular bentonite was then packed 
around the land surface to avoid preferential flow. The PVC pipes were capped to prevent direct 
rainfall and evaporation losses. The distance between Obs1 and Obs2 was approximately 300 m, 
while the distance between Obs2 and Obs3 was approximately 250 m. Obs1 and Obs2 were 
located in the low-lying area where the vadose zone consisted of fine to silty sand. In contrast, 
Obs3 was in the uplands, where both sand and sandstone lay within the unsaturated zone (Fig. 
1b). 

The hourly water table was monitored with a Keller DCX-22A data logger (Keller AG fiir 
Druckmesstechnik, Switzerland) that measured both the total pressure and barometric pressure, to 
account for the effect of barometric pressure. The barometric pressure transducers were 
positioned in a buffered thermal condition to exclude thermal effects (McLaughlin, 2011). A 
five-point moving average was used to smooth the groundwater level data (Loheide, 2008). 

A rain gauge (Onset Computer, USA) was installed to monitor hourly rainfall approximately 
400 m east of Obs3. Daily reference evapotranspiration (ETo) was calculated using an ETo 
calculator (Food and Agriculture Organization of the United Nations, Italy). The required 
meteorological data (i.e., net radiation, air temperature, air humidity, and wind speed) were 
collected from the closest meteorological station. 

The closest meteorological station is approximately 20 km southeast of the study area, but its 
monitoring data can represent the meteorological conditions of the study site because of its flat 
topography. Based on long-term data (2000-2017), the mean annual precipitation is 330 mm, with 
a significant annual variation indicated by a coefficient of variation (CV) of 0.3. Over 70% of the 
annual precipitation falls from June to September (Fig. 1c), and rainfall events with less than 5 mm 
are dominant, accounting for 72% of the total events (data not shown). The ET» calculated by the 
Penman-Monteith equation (Allen et al., 1998) is fairly stable at 1150 mm/a, with a CV as low as 
0.1. The monthly ETo ranges from 22 mm in January to 203 mm in May, with an average value of 
105 mm. The annual mean air temperature is approximately 7.0°C, with a maximum of 22.4°C in 
July and a minimum of —11.0°C in January (Fig. 1c). 
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Fig.1 (a) Location map showing the satellite image of the study sites (Obs1, Obs2 and Obs3); (b) the schematic 
hydrogeological cross section from measurements of the piezometers that are located at the bottom of wells; and 
(c) long-term average monthly precipitation, reference evapotranspiration (ETo), and air temperature 


2.2 Groundwater recharge estimation 


The WTF method is widely applied for estimating the groundwater recharge of shallow 
unconfined aquifers, assuming that the rise of the water table in unconfined aquifers results solely 
from rainfall infiltration (Healy and Cook, 2002; Varni et al., 2013). The method is best applied in 
areas with shallow water tables that demonstrate sharp rises in the water table over short time 
periods. A previous study (Yin et al., 2011b) indicated that the WTF method is applicable to this 
study site. Groundwater recharge (R) is calculated as follows: 

= am = eye : (1) 

dt At 

where Sy is the specific yield; d is the DWT (cm); h is the hydraulic head (cm); t is the time (h); 
Ah is defined as the difference between the peak of the rise and the low point of the extrapolated 
antecedent recession curve at the time of the peak to determine the total recharge; and At is the 
time span corresponding to Ah. 


2.3 ET, estimation 


Several methods are available for estimating ET, based on diurnal WTFs, including those based on 
the water balance principle (White, 1932; Loheide, 2008) and those based on statistical theory 
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(Soylu et al., 2012; Wang and Pozdniakov, 2014). The method developed by Loheide (2008) was 
adopted here because of its capability to estimate hourly groundwater consumption, rather than 
daily values commonly obtained from other methods. A method comparison conducted previously 
also indicated that the Loheide method yielded the most accurate estimation (Yin et al., 2013). The 
Loheide method is based on three assumptions: (1) the decline in groundwater levels is only due to 
plant root water uptake and surface evaporation; (2) ETg consumed by vegetation is negligible from 
midnight to 06:00; and (3) the rate of head change at the recovery source is equal to the overall rate 
of water table change at the observation location. The ET, can be calculated using the groundwater 
balance equation: 


Ah 

ET, ()=r()-S, re (2) 
where r(t) is the net inflow rate of groundwater (mm/h); Sy is the specific yield of sediments; and 
Ah/At is the change in DWT per unit time (mm/h). Because of the assumption that the rate of head 
change at the recovery source is equal to the overall rate of water table change at the observation 
location, the detrended water table depth (WTpr) is obtained by subtracting this trend from a 

portion of the observed water table records as follows: 
WT r(t) = WI(t)—m, xt—b,, (3) 
where WT(t) is the measured water table (mm); and mr and br are the slope (mm/h) and intercept 
of the trend (mm), respectively. A relationship can be determined to predict dWTpr/dt as a 


function, I (WTbpr), of the detrended DWT using periods of zero ET,. Subsequently, the net 
inflow rate of groundwater can be estimated using the following equation: 


r(t) =S, x[['(WT,,(¢)) +m, |. (4) 


As the fastest recovery of the water table occurred around 21:00 (LST) at Obs1 and 00:00 at 
Obs2 (Fig. 2a), the groundwater inflow was calculated using the recovery period from 21:00 to 
05:00 at Obs1 and 00:00 to 06:00 at Obs2, following the suggestion of Yin et al. (2013) that the 
fastest recovery of water table is an indicator of the suspension of ETg. 


2.4 Specific yield estimation 


It should be noted that Sy is the most crucial factor in determining the accuracy of recharge and 
ET (Nachabe, 2002; Loheide et al., 2005; Gribovszki, 2017). Sy is defined as the volume of water 
released under gravity from storage per unit cross-sectional area, per unit decline in the water 
table (Freeze and Cherry, 1979). When the water table is deep (>1 m), the application of a 
constant Sy is valid (Duke, 1972; Soylu et al., 2012). However, depth-dependent Sy should be 
considered in shallow water table conditions (Crosbie et al., 2005; Loheide, et al., 2005; Fan et 
al., 2014; Cheng et al., 2017) . At Obs1 and Obs2, with shallow water tables, two methods were 
employed to estimate Sy. The relationship between Sy and DWT can be described by Equation 5. 


nee gB 9, —, | (5) 
4 (1+ (ad)")™ 


where d is the DWT (cm); Os is the saturated soil moisture (cm?/cm°); 6, is the residual soil 
moisture (cm?/cm?); and a, n, and m are the empirical coefficients in the van Genuchten model 
(van Genuchten, 1980). Undisturbed soil samples were collected every 20 cm from the excavated 
profiles. Bulk density was estimated by the oven drying method, and particle-size distributions 
were determined by the sieving method for the sand fractions and by a soil hydrometer for the silt 
and clay fractions, respectively. 

Sy was also determined by the water table response to rainfall events, as the ratio of infiltrated 
rainfall to recorded rise in water table level. Rainfall events were only considered when soil 
moisture was likely to be replenished by previous rainfall. The interception loss by the plant 
canopy was ignored because the grasses were small and sparsely distributed at the study sites. The 
Sy of sandstone at Obs3 was assumed to be constant because the water table was deeper than 3 m. 
We assigned the value of 0.08 to S, based on a previous study (Yin et al., 2011b), which 
determined it through a pumping test that used the type-curve method. 
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3 Results and discussion 


3.1 Variation of water table 


The monitored DWT values at Obs1 and Obs2 from December 2017 to March 2018 were 
excluded because the soil was frozen. The patterns of water table variations were strongly related 
to DWT. At Obs1, Obs2, and Obs3, DWT varied between 0 and 70, 0 and 100, and 330 and 370 
cm, respectively (Fig. 2a), with mean values of 20, 60, and 350 cm, respectively. Hourly 
groundwater levels at Obs1 responded to rainfall rapidly (usually within 1 h, data not shown), and 
levels were the highest among the three observation sites. The water table at Obs2 had similar but 
somewhat slower responses to rainfall (within 1-2 h) and was lower compared to that at Obs1. 
For example, the water table raised by 10 cm following a rainfall of 11.0 mm/d on 8 July 2017, at 
Obs2, while it rose by 27 cm at Obs1. The water table at Obs3 displayed a smooth and seasonal 
variation, and only reacted to heavy rainfall at the end of July and October, with a delay of 
approximately four days (Fig. 2a). 

Diurnal fluctuations were observed for the entire observation period at Obs1, except during 
rainy periods (Fig. 2a). At Obs2, diurnal signals were obvious when the DWT was less than 75 
cm, and diminished when the water table was lower, a behavior previously observed by others 
(Mould et al., 2010; Nippert et al., 2010). The diurnal magnitude was much higher at Obs1 than at 
Obs2, which had a deeper water table. For instance, the mean magnitude was 4.4 cm during 5-14 
June 2018 at Obs1, three times higher than that at Obs2 (Fig. 2a). It was also observed that the 
water table began to decline approximately 2 h earlier at Obs1 than at Obs2 (Fig. 2b). In contrast, 
the water table at Obs3 did not show any diurnal fluctuations (Fig. 2a). The response to rainfall 
and the magnitude of diurnal fluctuations are attributed to the thickness of the unsaturated zones, 
which determines whether plants can extract groundwater. 

The seasonal water table regime was controlled by meteorological conditions (rainfall and ET»). 
In the study area, a hydrological year can be divided into a dry season (April—mid-July) and a wet 
season (the remaining period) (Yin et al., 2015). The DWT at Obs1 and Obs2 also reflected the 
characteristics of dry and wet seasons (Fig. 2c). In the dry season, the water table generally 
declined in response to lower rainfall and high ETo, while in the wet season, the water table rose 
because of the increase in rainfall and decline in ET». 
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Fig. 2 Water table regimes at three observation sites, Obs1, Obs2, and Obs3 (a and b), daily rainfall and 
potential evapotranspiration (ETo) (c). DWT, depth to water table. 


3.2 Specific yield (S,) 


Soil texture analysis indicated that the soil was quite homogenous at Obs1 and Obs2; therefore, a 
uniform texture was assumed to estimate S, using Equation 5. Sy estimated using the water table 
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response to rainfall and by soil texture, decreased with DWT (Fig. 3). The sample data and fitted 
exponential function (R*=0.75, P<0.01) are shown in Figure 3. For silty sand soils, such as those 
at our sites, Sy was found to increase with DWT and also fitted an exponential pattern, reaching a 
relatively stable value within 1 m (Crosbie et al., 2005; Fan et al., 2014; Hill and Durchholz, 
2015; Cheng et al., 2017). The relationship indicated that S, varied significantly, with the 
difference of an order of magnitude in extreme cases, suggesting that ignoring Sy variations will 
result in substantial interpretation errors at sites with a shallow water table. 

However, it should be noted that Sy from the water table response method was lower than that 
obtained using the soil texture method (Fig. 3). Differences have been widely observed among Sy 
obtained from commonly used methods, that is, soil texture, drainage, and water table response 
methods (Fan et al., 2014; Gribovszki, 2017). Method comparisons indicate that S, values 
obtained using the water table response method are suitable for analyzing groundwater recharge, 
as the same temporal and spatial scales are used (Crosbie et al., 2005). Therefore, S, from the 
water table response method was applied to estimate the groundwater recharge, whereas Sy from 
the soil texture method was used to calculate ETs. 
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Fig. 3 Increases in specific yield (Sy) with DWT using the water table response method and the soil texture 
method 


Soil water moisture before rainfall is a vital factor when using the water table response method 
to estimate Sy. Rainfall recharges soil water if the water content of the unsaturated zone is low, 
leading to an overestimation of Sy. Therefore, such rainfall events were only considered when soil 
moisture in the unsaturated zone was likely to be replenished by previous rainfall within one 
week of S, estimation (Fan et al., 2014). However, this can also lead to significant errors when the 
water table is deeper, similarly to the Sy with DWT of 92 cm in Figure 3. In addition, intense 
rainfall can cause air from the unsaturated zone to infiltrate into the saturated zone, leading to a 
rapid rise in the water table. This process is called the Lisse effect (Crosbie et al., 2005; Miyazaki 
et al., 2012), which often causes the underestimation of Sy. However, in the study area, the Lisse 
effect was considered to be minimal, owing to the coastal soil texture. 


3.3 Groundwater recharge 


For Obs1 and Obs2, groundwater recharge was estimated for 26 rainfall events during June 2017 
to June 2018. Groundwater recharge generally increased with the daily rainfall intensity (Fig. 4a). 
The rate of increase slowed when rainfall exceeded 20 mm/d and even decreased when rainfall 
was higher than 30 mm/d at Obs2. This phenomenon was attributed to heavy rainfall exceeding 
the infiltration capacity of soil (Horton, 1933; Beven, 2004; Liu et al., 2011), or to the elevation 
of the water table to the land surface, leaving no space for infiltration (Sophocleous, 2002; Hoang 
et al., 2017), with rainfall becoming surface runoff. 

At the event scale, groundwater recharge was also affected by hourly rainfall intensity. For 
example, the rainfall amount was approximately 28.6 and 33.2 mm on 24 and 26 July 2017, 
respectively, however, the groundwater recharge on 26 July was much lower (13.4 mm) than that 
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on 24 July (23.9 mm) at Ob2 (Fig. 4b). Although rainfall was heavier on 26 July, it mostly fell 
within 1 h (approximately 21.7 mm/h). In comparison to the prolonged rainfall on 24 July, the 
concentrated rainfall reduced groundwater recharge, as rainfall mainly became surface runoff. 
The relationship between rainfall intensity and groundwater recharge indicates that groundwater 
resources will decrease under the projected climate change, in which heavy rainfall events will 
increase with an overall reduction in rainfall frequency (Allen and Ingram, 2002). 

Antecedent rainfall also affected groundwater recharge in some cases. For instance, 
groundwater recharge was approximately 6.9 mm/d when the rainfall amount was 14.9 mm on 18 
June 2017, while it reached 12.8 mm on 21 June 2018 when the rainfall was slightly lower (13.1 
mm/d). It was found that the five-day cumulative rainfall was 1.5 mm for the first case and much 
lower than 10.2 mm for the latter event. The antecedent rainfall increased soil moisture, and 
subsequent rainfall replenished soil moisture and produced more recharge than the case in which 
the soil was dry and most infiltrated rainfall was retained in vadose zones (Zhang and Schilling, 
2006). 

In addition, groundwater was also determined by DWT as groundwater recharge at Obs2 was 
higher than that at Obs1, with a shallower water table in most cases (Fig. 4a). However, the 
relationship between groundwater recharge and DWT was nonlinear. For example, the total 
groundwater recharge during the study period was the highest at Obs2, corresponding to a mean 
DWT of 72 cm in the rainy period, compared to that at the shallowest site (mean DWT, 30 cm at 
Obs1) or at the deepest site (mean DWT, 360 cm at Obs3). It seems that both deeper and 
shallower water tables restrict groundwater recharge. In areas where the water table is close to the 
land surface, such as Obs1, rainfall mainly becomes surface runoff due to the saturation of the 
ground, instead of infiltration to groundwater (Gillham, 1984). At sites with thick unsaturated 
zones, such as Obs3, the majority of rainfall infiltration is consumed by evapotranspiration in the 
uppermost 2 m of soil in semi-arid regions (Cook et al., 1989), and only a small portion reaches 
groundwater. 

At the monthly scale, the temporal variations of groundwater recharge at Obs1 and Obs2 were 
similar and generally determined by the amount of rainfall (Fig. 4c). The greatest recharge 
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Fig. 4 Relationships between event-based groundwater recharge and rainfalls at Obs1 and Obs2 (a), water table 
fluctuations in response to rainfall at 00:00 on 24 and 26 July 2017 (b), monthly groundwater recharge and 
rainfall (c), and their fitted relationships (d) 
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occurred in July 2017, corresponding to the heaviest monthly rainfall of 94.2 mm, while the 
lowest recharge occurred in September 2017, when the monthly rainfall was as low as 16.6 mm. A 
linear equation with a positive slope can be used to describe the relationship between the monthly 
groundwater recharge and rainfall (Fig. 4d). Groundwater recharge was usually higher at Obs2 at 
a monthly scale than that at Obs1 (Fig. 4c), as the mean DWT for Obs1 and Obs2 was 21.0 and 
61.0 cm, respectively. During the study period, total groundwater recharge was approximately 
207.0 mm, accounting for 52% of the rainfall at Obs1, while groundwater recharge was 
approximately 250.0 and 21.0 mm at Obs2 and Obs3, respectively, accounting for 63% and 5% of 
rainfall, respectively. 

Overall, the recharge rate of Obs2 was higher than those of Obs1 and Obs3, indicating that both 
shallower and deeper water tables can limit rainfall recharge. The soil texture of the unsaturated 
zone in Obs3 was also a major factor that decreased the recharge rate. Compared with Obs1 and 
Obs2, similar recharge rates (53% and 62%) were also found in the study by Crosbie et al. (2005) 
and Fan et al. (2014) with similar soil textures and plants. Although the mean DWT (120.0 and 
200.0 cm) of their study areas was greater, the annual rainfall was almost four times higher than 
that at our study site. Therefore, a greater DWT does not decrease the rainfall infiltration in case 
of a high water content in the unsaturated zones. 


3.4 Groundwater evapotranspiration 


As no diurnal fluctuations were observed at Obs3 (Fig. 2a), hourly ET, (excluding rainy days), 
was estimated for Obs1 and Obs2 only. The period from November 2017 to March 2018 was 
excluded because the ET, values were negligible due to plant dormancy (Yin et al., 2014). 
Therefore, ET, was calculated for 129 d for Obs1 and 133 d for Obs2. It was observed that ET, 
clearly varied diurnally. The mean ET, for each hour shows that ET, usually started increasing 
from 06:00 to a maximum of approximately 0.54 mm/h at 13:00 to 14:00, then decreased back to 
zero at 21:00 at Obs1 (Fig. 5). From 22:00 to 05:00 the next day, ET, was close to zero, indicating 
that there was no groundwater uptake during the night. Hourly ET, variations at Obs2 were 
similar to those at Obs1 (R’=0.91, P<0.01), except that the peak values were lower 
(approximately 0.29 mm/h) and the duration of the peak was longer (approximately 4 h versus 2 h 
at Obs1) (Fig. 5). During daytime, the standard deviation (SD) of ET, increased with time, 
reaching a peak value at 13:00, and then decreased. For example, SD increased from 0.05 mm/h 
at 06:00 to 0.27 mm/h at 13:00, and decreased nearly to 0.00 mm/h at 21:00 at Obs1. 
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Fig.5 Mean hourly ETg (groundwater evapotranspiration) at Obs1 (a) and Obs2 (b) 


Daily ET, patterns suggested that it started to increase around mid-April when plants began to 
grow and became nearly zero in mid-October when they entered winter dormancy (Fig. 6), with 
good agreement on daily sap flow variations in the Mu Us Sandy Land of China (Yin et al., 2014). 
At a daily scale, ET, at Obs1 ranged from 0.3 to 8.7 mm/d, with a mean value of 4.3 mm/d, while 
at Obs2 it varied between 0.1 and 7.8 mm/d, with a mean value of 2.5 mm/d. It should be noted 
that ET, at Obs1 and Obs2 was in line with other grass-covered areas, where it is reported to be in 
the range of 1.0—4.0 mm/d (Nachabe et al., 2005; Schilling, 2007; Lautz, 2008; Gribovszki et al., 
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2017). Meteorological conditions were the main factors affecting the fluctuations in daily ET,. 
The correlation analysis revealed that net radiation and vapor pressure deficit and relative 
humidity were positively and negatively related to ET,, respectively (Fig. 7a, b and c). However, 
there was no correlation between ET, and wind speed (Fig. 7d). The most influential 
meteorological factors on ET, were net radiation, relative humidity, vapor pressure deficit, and 
wind speed. 

Although hourly and daily ET, at Obs1 and Obs2 exhibited similar temporal variations, the 
magnitudes were substantially different (Figs. 5 and 6). It was hypothesized that DWT, the main 
difference between the two sites, caused this difference. To explore the impact of DWT on ET,, 
we introduced the ratio of ET, to ETo, as meteorological conditions also affect ET, in shallow 
water table conditions, as shown in Figure 7 and other studies (Butler et al., 2007; Carlson Mazur 
et al., 2014; Shanafield et al., 2017). ET,/ETo gradually decreased with DWT and an exponential 
equation was used to fit this decay (Fig. 8a), which is consistent with field observations (Nichols, 
1994; Cooper et al., 2006; Elmore et al., 2006; Pritchett and Manning, 2012; Ma et al., 2013; 
Shen et al., 2015; Liu et al., 2017) and numerical modeling (Shah et al., 2007; Soylu et al., 2012). 
It is predicted that, based on an extrapolation using the fitted equation, ET,/ETo will drop to 0.005 
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Fig. 7 Relationships of daily ET, with net radiation (a), vapor pressure deficit (b), relative humidity (c), and 
wind speed (d) 
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when the water table decreases to 3.0 m. The depth corresponding to an ET,/ETo value below 
0.005 is defined as the extinction depth (Shah et al., 2007), which is a key parameter for 
quantifying regional evapotranspiration. The extinction depth of the study site was approximately 
3.0 m, explained by the absence of diumal fluctuations at Obs3, where the mean DWT was 
approximately 3.5 m. 

The monthly ET, was obtained by summing daily values. On a monthly scale, ET, ranged from 
125.0 mm in June 2018 to 4.0 mm in April 2018, with a mean value of 67.0 mm at Obs1, while 
ET, ranged from 94.0 mm in May 2018 to 8 mm in October 2017, with a mean value of 42.0 mm 
at Obs2. The coefficient of determination (R?) between ET, and ETo at Obs1 was 0.94, which was 
higher than that at Obs2 (R?=0.29), indicating that ET, at Obs1 was mainly determined by 
meteorological conditions. ET, at Obs2 was also controlled by DWT, as indicated by the lower 
values in July 2017 and June 2018, when the water table was at its lowest (Figs. 2a and 8b), 
although ETo in those times was higher than that during the neighboring months (Fig. 8b). The 
monthly variations in ET, suggested that ET, was determined by the DWT and meteorological 
conditions. 

Over the hydrological year studied, the groundwater consumption at Obs1 was approximately 
539.0 mm, which was 1.6 times the annual rainfall. At Obs2 it was approximately 336.0 mm, 
equivalent to 0.9 times the annual rainfall, indicating strong groundwater uptake by ET, in 
shallow water table areas. In shallow water table environments, ET, can significantly exceed 
rainfall, as shown in this and previous studies. For example, in a playa in arid central Australia, 
annual ET; loss is as high as annual rainfall (Shanafield et al., 2015). In the riparian zones of the 
Middle Rio Grande River, the estimated ET, is even higher, approximately four times the annual 
rainfall (Soylu et al., 2012). 

Additionally, the variation patterns of rainfall recharge and ET, were only analyzed for one 
hydrological year because of limited monitoring data. In the future, the interannual variations of 
rainfall recharge and ET, and their determining factors should be explored. 
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Fig. 8 Relationship between DWT and ET,/ETo (a) and monthly ETo and ETs variations at Obs1 and Obs2 (b) 


4 Conclusions 


In this study, WTFs in three observation wells along a cross-sectional area were observed over a 
hydrological year, in a semi-arid lake wetland in the Ordos Plateau, China. The data were 
analyzed to estimate groundwater recharge and evapotranspiration using a depth-dependent 
specific yield. A rapid water table response to rainfall was observed for two shallow wells (i.e., 
Obs1 and Obs2), at which diurnal WTFs were also observed over most of the study period. This 
indicated that the most influential factors of variations in water table level are rainfall and ET,. 
However, at the deeper water table site (Obs3), the groundwater table only reacted to extremely 
heavy rainfall and showed no diurnal variations. 

Over the study period, the annual groundwater recharge was 207.0, 250.0, and 21.0 mm, 
accounting for 52%, 63%, and 5% of rainfall at Obs1, Obs2, and Obs3, respectively. It was found 
that rainfall intensity is one of the main determining factors for rainfall infiltration, as 
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groundwater recharge increases with rainfall intensity at daily and monthly scales, and the rate of 
increase of daily groundwater recharge slowed down when rainfall exceeded 20-30 mm/d, due to 
surface runoff. In addition, DWT was also a major factor influencing groundwater recharge, 
according to the recharge rate of the three observation wells in this study. 

Hourly ET, showed that groundwater consumption occurred mainly during the daytime, and 
nighttime ET, was negligible. The daily mean ET, at Obs1 and Obs2 were 4.3 and 2.5 mm, 
respectively, over the entire hydrological year. The variations in ET, were mainly attributed to the 
DWT and net radiation in the study area. 
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